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Abstract 

It is usually believed that the spin-fluctuation 
mechanism for high-temperature superconduc- 
tivity results in d-wave pairing, and that it is de- 
structive for the conventional phonon-mediated 
pairing. We show that in bilayer materials, due 
to nearly perfect antiferromagnetic spin correla- 
tions between the planes, the stronger instabil- 
ity is with respect to a superconducting state 
whose order parameters in the even and odd 
plane-bands have opposite signs, while having 
both two-dimensional s-symmetry. The interac- 
tion of electrons with Raman- (infrared-) active 

phonons enhances (suppresses) the instability. 
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the symmetry of the pairing state [|[]. In- ers. In this aspect our mechanism is very 
timately related to this, is the question of different from the interlayer j^azr-tunneUing 
whether the superconductivity is due to an- (IPT) mechanism discussed by Anderson et 
tiferromagnetic spin fluctuations (see e.g. al. [Q. Nevertheless, similar to the IPT 
Monthoux and Pines (MP) Ref. [Q, and also model, any attractive interaction between 
Refs. to electron-phonon (EP) interac- electrons in the same band, such as the one 
tion enhanced by inter-layer pair tunnelling mediated by even (Raman-active) phonons, 
or to neither of two. In this discus- enhances T^. This is opposite to the previ- 
sion, it is indirectly assumed that the antifer- ously considered single-layer AFSF models in 
romagnetic spin-fluctuation (AFSF) mecha- which such interactions are mutually destruc- 
nism necessarily leads to d-wave pairing, and five. 

that the AFSF and EP mechanisms cannot Support for such an enhancement mech- 

coexist. anism may be found in the experimental 

In this Letter we point out that, whereas fact (e.g. Refs. PJTO[|) that some mem- 

the AFSF mechanism leads to rf-pairing for bers of the cuprate family (Nd2_xCea.Cu04, 

one layer, it may lead to (two-dimensional) HgBa2Cu04) behave as conventional s-wave 

s-symmetry for a bilayer. The condition for EP superconductors. MP AFSF theory, on 

that is existence of strong antiferromagnetic the other hand, would have to imply princi- 

correlations between the two layers in a bi- pally different mechanisms for this compound 

layer, as found experimentally in YBa2Cu307 and for those with high Tc's. Another ex- 

PJ^Q. We find that, for a given coupling perimental fact which suggests a constructive 

strength, Tc(s) is about twice as high as interplay between phonon- and non-phonon 

Tc{d) thus making it easier to achieve the mechanisms is that in YBa2Cu307, the iso- 

observed values of Tc ~ lOOK. Essential for tope effect increases smoothly when the su- 

the positive influence of layer-doubling is the perconductivity is suppressed |1T[ . Finally, 



sm^f/e-particle tunneling which splits the one- the most impressive argument is that in all 
electron plane-bands into even and odd with high-T^ materials is anticorrelated mth the 
respect to the mirror-plane between the lay- in-plane antiferromagnetic correlation length 



^. In particular, in YBa2Cu307.o, ^ is about regards the interplane hopping inside the unit 

one lattice parameter, which would make the cell, t± (k), we assume that it is sufficiently 

single-layer AFSF superconductivity virtu- large to set even and odd symmetry of the 

ally inoperative. To the contrary, as we shall two-dimensional bands, but we neglect, for 

argue below, the proposed bilayer model is simplicity, in the following numerical calcula- 

barely sensitive to the in-plane AF correla- tions the even-odd splitting e_ (k) — e+ (k) = 

tion length ^ at all. 2i:^(k). 

In the following we shall assume a con- The generalization of the MP AFSF 

ventional picture in the sense that the one- model to two-bands is straightforward; one 

electron tunneling between the planes is al- has only to take into account that the effec- 

lowed both in the normal and in the super- tive vertex for scattering of an electron from 

conducting states. This is in contradiction band i to band j by a spin-fluctuation de- 

with the IPT scenario ||^, but in agreement pends on i,j, while the spectrum of the fiuc- 

with some photoemission experiments [|12|. tuations x is the same as in MP. Then, Eqs. 

In this case, the single-particle eigenstates for (6-8) of MP become 



a bilayer are the even |+,k) and odd I— ,k) \^ /i ■ n t^v^v^t/ /i • • n 
combinations of the individual plane states 

X Gki (q, lu), 

and k is the 2D Bloch-vector. The properties 



(k, iun) = [iujn - e (k) + /i] 6ij - (k, iu^) 
of the even and odd bands are discussed in de- 



^ij (k, iuJn) = ^'^^ij ~ 



tail in Ref. ||l4[, but, for the purpose of com- ust 



liOn — HjJ„ 



parison with the MP model, we use the same ^ ^' ' ' 
band model as they did. We neglect com- where S and $ are respectively the normal 
pletely the kz dispersion due to small inter- and anomalous self-energies, G is the single- 
cell c-hopping, which can lead to interesting particle Green function, e the bare electron 



effects (see, e.g., Ref. [0) but which are how- energy, and I]qm denotes the average over 
ever beyound the scope of this Letter. Ac- the Brillouin zone plus the sum over the Mat- 
cordingly, in the following the term "bands" subara frequencies. V is AFSF pairing inter- 
always means "two-dimensional bands". As action, determined by the exchange interac- 



tion of electrons with the AFSF's, Vij, kl = S_ (k, icUn) = Tg"^ ^ x (k - q, icUn - itOm) G+ (q, icUm) 

qm 

JdRdR'j:ap^5MJir-R)a^^\jp) G-^{k,tuj^) =^cu,-e(k) +/x-S+(k,2^„) 

X m-^'){kl\Ar-R)a,s\l6) where J is ^j,, = -T/ ^ X (k - q, - z^m) (2) 

exchange interaction and x = ('S(R.)S(R')) '^"^ 

X (-q, -iuJm) ^- (q, iuJm) G- (q, zu;„^) 

is spin-spin correlation function. For a bi- 

layer, one can let R be a two-dimensional vec- and the same with + and - subscripts in- 

tor and introduce x = - ^')Iuv, where terchanged, and g and x are the same as in 

u,v = 1,2 label layers, and / accounts for in- MP. 

terplane correlations, if any. Then the func- For reasons of symmetry, the solution of 
tion X is the same as in MP. these equations must have the form G+ = 
The key to the our bilayer AFSF model $+ = ±$_. For the upper choice of the 
is the experimental fact that the spin flue- sign, the Eqs. (D reduce precisely to the orig- 
tuations in the bilayer of YBa2Cu307_:j; are inal MP pairing state. For the lower choice, 
always antiferromagnetically correlated be- the Eqs. (|) again reduce to the one-plane 
tween the planes ||,?,||. Even fully oxy- case, but now the interaction m the equa- 
genated samples, where the in-plane corre- i^on for $ is effectively attractive. In other 
lation length is already of the order of the words, now the order parameter has the op- 
lattice parameter, show nearly perfect inter- posite sign in the two bands, and therefore 
layer correlation g]. The exchange potential the last Eq. in (|) can be rewritten in terms 
set up by such a spin-fluctuation is there- of |$|, and with plus instead of minus on the 
fore odd with respect to the mid-layer mir- right-hand side. 

ror plane and, correspondingly, couples ex- The concept of a superconducting state 

clusively even and odd electron states (but where two distinctive bands had the order 

neither odd to odd, nor even to even). In parameters of the opposite signs was first dis- 

other words, = -J„„ = -1. In this case, cussed in 1973 in connection with semimet- 

after summation over u,v in the expression als More recently, in a two-layer Hub- 



for Vij-fc/ and defining the appropriate cou- bard model, such a solution was found by 
pling constant g, Eqs. become: Bulut et al |T3 (which they labeled as "4" 



state) and in the conventional superconduc- more is lost by making an s-state than by 

tivity theory ||18| , where it appears in case of making a d-state (which has been found nu- 

strongly anisotropic electron-phonon and/or merically by MP), because the latter loss is 

Coulomb interaction, or because of a strong the difference between the small-g loss and 

interband scattering by magnetic impurities, the large-g gain, while in an s-state one loses 

In all cases, order parameter has s-symmetry over the whole Fermi surface . 



inside each band and changes sign between Now, coming to the bilayer case, we ob- 
the bands. serve that there is no conflict between the 
From Eqs. @ it is quite plausible that small and the large g's any more. The 
such an instability is stronger than the dr^2_y2 AFSF interaction spans two different sheets 
one, and will occur at a higher Tc. Below we of the Fermi surface, which always have or- 
shall prove this numerically, but before go- der parameters of the opposite signs. Thus 
ing to numerical results, it is instructive to the AFSF interaction is as attractive for s- 
get a conceptual understanding about these pairing in a bilayer as it is repulsive in a sin- 
two different solutions. The physical reason gle plane, and consequently more attractive 
for having rf-symmetry in the one-plane case, than (i-pairing in a single plane. Of course, 
is that the AFSF interaction makes pair- the resulting s-state is likely to be highly 
ing energetically favorable only when it cou- anisotropic, to take better advantage of the 
pies parts of the Fermi surface which have large x(ct) at q ^ Q. This is similar to An- 
opposite signs of the order parameter [|l9l. derson's model 0. To demonstrate this ef- 
In Y123 the AFSF interaction is peaked at feet numerically, we have solved Eqs.|^ with 
Q=(7r/a, 7r/a). The shape of the Fermi sur- the parameters from MP paper, and using 
face is such that the condition is satisfied only the same numerical technique. As expected, 
for dj.2_y2 symmetry. On the other hand, the the maximal eigenvalue of the last Eq.0 is 
small-q interaction couples parts of the Fermi larger than that for the MP rf-pairing (about 
surface where the order parameter has the 1.5 compared to 1). Fig.|l] shows the plot of 
same sign. This makes pairing unfavorable. as a function of the interaction constant 
Since x(q ~ Q) ^ x(q ~ 0), nevertheless, g for both cases. To test the numerics, we 



have also solved the original MP equations frequency-dependent superconducting gap is 
and obtained the similar results as MP. related to $ as 

From Fig.|I|, one immediately observes $(p,^w„) $(p,^cj„) 

A{p,lUJn) 



that the value (0.69 eV) of the coupling con- ' ^(P' ^^") ^ " /mS(p, ia;„)/u;„ 

stant g which yields ~ 90K for two planes From Fig.^we observe that the absolute value 

and s-symmetry, is much smaller than the of A behaves similarly in both cases, having 

corresponding value (1.24 eV) for one plane a minimum along (11) directions and a max- 

and (i-symmetry. Actually, Tc(s) ~ 2Tc((i) imum along (10) directions. Furthermore, 

for g up to about 1 eV. At stronger cou- |A|'s in both cases differ by less than 10% on 

plings, Tc{d) saturates faster than Tc(s) due two-thirds of the whole Fermi surface, thus 

to stronger effect of mass renormalization. making it extremely difficult to distinguish 

Similarly, as we shall see below, the ratio of between the two in an experiment which does 

the maximal gap to tends to be larger for not probe the relative phases of A. In other 

the one-plane model, for the same Tc. words, the order parameter, formally having 

One can also obtain the self-consistent so- s— symmetry, is still strongly anisotropic, but 

lution for $ at T <^ Tf.. To do that, one has nodeless. 

to include higher-order terms (see, e.g., Ref. Now we shall briefly discuss some exper- 

0]). In this case one of the Green functions imental consequences of the bilayer AFSF 

in the Eq.^ for $ should be replaced by: superconductivity model. It turns out that 

~ many difficulties associated with the original 

G_ {k,lUJn) = tUJn - (Ck - ;u) - S_(k,zt^„), (3) 

AFSF superconductivity model disappear in 



the present version. 



— $-(k, iuJn)G-{—'k, — ico'„)<l>_(— q, — ic<j„)]. 

(1) In-plane vs. perpendicular-to-the- 

This new set of equations can be solved iter- planes Josephson tunneling. Recent searches 



atively, starting with G = G (which is cor- for the rf-pairing in YBa2Cu307 (Refs. pT|j22 

rect to first order in $). The actual solution and others) still do not give a definite an- 

for T = Tc/2, shown in Fig.Q, was achieved swer. Experiments probing the angular de- 

by making two iterations of EqsJ^. The pendence of the order parameter, as well 



the existence of the so-called "paramagnetic completely decoupled from the planes, but 

Meissner effect" |2^, indicate the existence of this is inconsistent with the strong effect on 

order parameters of opposite signs; Many ex- the Ba Aig phonon mode of the onset of su- 

periments were interpreted in terms of (ii,2-j;2, perconductivity). This finds natural expla- 

but such an interpretations can be questioned nation in our AFSF model: The chain im- 

because of presence of the chain band DTSf. purity potential is even with respect to the 

On the other hand, the existence of the finite mid-plane reflection and does therefore not 

tunneling current perpendicular to the planes produce scattering between the even and odd 

^ is incompatible with (i-pairing, but com- bands. 



patible with our model once the simplifying (4) The MP model has difficulty in ex- 
assumption e_ (k) = £+ (k) is removed p4[. plaining the continuous change of the isotope 

(2) The original MP model is very sensi- effect with oxygen content, as well as in rec- 
tive to the in-plane correlation length ^ : Ac- onciling the apparent phonon s— wave super- 
cording to Ref . [0] , Tc drops from 90K to 35K conductivity in Nd-cuprate with the assumed 
when ^ is reduced from 2.3a to a. Experimen- AFSF superconductivity in YBa2Cu307. The 
tally, in fully oxygenated samples (Oe.g-T.o) C, basics of this conflict is as follows: The AFSF 
is small and it still decreases closer to the interaction is pairing when a pair changes the 
Oy-composition where it becomes less than sign of its order parameter upon scattering; 
a. Tc is, however, not sensitive to the oxy- this is the case MP model for q ^ (1, l)7r/a. 
gen content in this regime. In our model and It is Repairing if there is no sign change. The 
contrary to MP, is hardly sensitive to the opposite is true for the electron-phonon in- 
sharpness of x{(l)> which predominantly in- teraction. Obviously, the only way to rec- 
fiuences the gap anisotropy. oncile the MP model with the known facts 

(3) One of the arguments in favor of the about the role of phonons in superconduc- 
AFSF mechanism has been the strong Tc sup- tivity is to assume that the electron-phonon 
pression upon Zn doping. However, it has re- interaction, contrary to the AFSF interac- 
mained unclear why chain disorder hardly af- tion, is strong for q — and weak for 
fects Tc. (One could argue that the chains are large g's; this is opposite to common wis- 



dom. In our model the corresponding as- 
sumption is much less painful: One has to 
assume that the even phonons, like for in- 
stance Aig Raman-active phonons, interact 
with electrons stronger than the odd ones. 
This seems quite plausible and some indirect 
arguments can be given in support of this as- 
sumption (e.g., the even phonons strongly in- 
fluence the extended van Hove singularities 
\U . 

In conclusion, we have shown that from 
the observed strong antiferromagnetic cor- 
relations between the Cu-0 planes in bi- 
layer materials the strongest superconduct- 
ing instability due to antiferromagnetic spin 
fluctuations appears in the anisotropic s- 
channel, but so that the order parameters 
in the bonding (symmetric) and antibonding 
(antisymmetric) bands have opposite signs. 
This helps to reconcile the magnetic-induced 
superconductivity model with many exper- 
iments which previously seemed to contra- 
dict the magnetic scenario. In particular, the 
interrelation between the doping dependen- 
cies of the magnetic and the superconducting 
properties can be much easier understood. 

We want to thank O. Gunnarsson for 
many helpful discussions, especially regard- 



ing the relation of this model to the two-band 
Hubbard model. 
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FIGURES 

FIG. 1. Critical temperature as a func- 
tion of the coupling strength g for a single 
plane ((i3,2_y2)-symmetry) and for two planes 
(p^-symmetry). Parameters are as in Ref. [||. 

FIG. 2. Superconducting gap A for the low- 
est Matsubara frequency ujn = vrfc^T in the 2D 
Brillouin zone at T = 0.5Tc for a single plane (a), 
and for two planes (b). The Fermi-surface con- 
tour shown at the base of the plots corresponds 
to Ef = -0.37 eV (0.25 holes/plane). 
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